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The baby Skyrme model is a well-known nonlinear field theory supporting topo- 
logical solitons in two space dimensions. In the limit where the term quadratic in 
derivatives (the "sigma model term") vanishes some additional structure emerges. 
The resulting ("extreme" or "restricted" or "BPS") baby Skyrme model has exact 
soliton solutions saturating a BPS bound which exists for this restricted model. Fur- 
ther, the restricted model has infinitely many symmetries and infinitely many con- 
servation laws. Here we consider the gauged version of the restricted baby Skyrme 
model with gauge group U(l) and the usual Maxwell term for the gauge field. We 
find that, again, there exists a BPS bound and BPS solutions saturating this bound. 
^IIh ' We further find that the whole problem is essentially determined by a new kind 

of superpotential equation. The BPS bound and the corresponding BPS solitons 
J> , only may exist for potentials such that the superpotential equation has a solution 
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which exists globally, i.e., on the whole target space. We also calculate soliton solu- 
tions both exactly and numerically, completely confirming our qualitative analytical 
results. 
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PACS numbers: 11.30.Pb, 11.27.+d 



^ 1 I. INTRODUCTION 

It is a notoriously difficult problem to derive the basic physical properties of hadrons 
and nuclei from the underlying fundamental theory of strong interactions, i.e., quantum 
chromodynamics (QCD). One important approach to resolve this problem consists in the 
introduction of low-energy effective field theories, where the fundamental fields of QCD, 
which are the relevant degrees of freedom at high energies, are replaced by some other 
degrees of freedom which are supposed to describe strong interaction physics at low energies 
more practically or more concisely. The Skyrme model JlJ is one well-known proposal for 
such an effective field theory, where the basic degrees of freedom consist of a triplet of chiral 
fields (Skyrme fields) which are related to the Goldstone bosons of chiral symmetry breaking 
in the underlying QCD. The nucleons and nuclei, on the other hand, emerge as collective 
excitations of the basic fields, i.e., topological soliton solutions (called skyrmions) which the 
theory supports. Further, the topological degree of the solitons is identified with the baryon 
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number. Quantization of the spin and isospin degrees of freedom for the simplest solitons 
leads to predictions for the physical properties of nucleons and light nuclei which are in 
reasonable agreement with experiment, see, e.g., j^J, jsj, 0] (for further, more recent results 
see, e.g., jsj). These early successes have lead to an extensive numerical study of soliton 
solutions, both for the standard Skyrme model (consisting of a non-linear sigma model term 
quadratic in first derivatives and a so-called Skyrme term which is quartic in first derivatives) 
and for the standard Skyrme model with an additional potential term ("pion mass term"). 
As a result, solitons in the standard Skyrme model for topological degree up to about 20 
are now fairly well understood 0], 0], 8j. 

Although the Skyrme model lagrangian cannot be rigorously derived from QCD, the 
symmetries and anomalies of QCD provide sufficient information to uniquely determine the 
coupling of the Skyrme field to electromagnetism 9]. Concretely, the gauge field couples 
both to the topological (baryon) current and to the third component of the isospin current. 
A first attempt to determine the electromagnetic contributions to the proton and neutron 
masses within the Skyrme model perturbatively in the electromagnetic field and within 
the Couloumbic approximation has been done in [loj|. Unsurprisingly, the electrostatic 
energy of the (charged) proton is larger than that of the (electrically neutral) neutron. A 
correct calculation of the proton-neutron mass difference most likely requires some amount 
of explicit breaking of the isospin symmetry. A mainly numerical analysis of the standard 
Skyrme model with the electromagnetic field coupled to the third component of the isospin 
has been performed in ll|, and some further results have been presented in 12]. 

To the best of our knowledge, the analysis of the fully nonlinear coupled Skyrme-Maxwell 
system has not yet been carried to the point where reliable descriptions of the electromag- 
netic properties of nucleons and light nuclei would be possible, despite the obvious physical 
interest of such a description, most likely due to the tremendous difficulty of this task. It 
is at this point where the study of lower- dimensional models becomes useful or even indis- 
pensable. The investigation of lower dimensional field theories is a well-established research 
field, both in order to understand and solve difficult nonlinear problems, and because of 
the growing appearance of interesting physical realizations of planar physical systems. The 
simplifications which may be achieved in lower dimensions help to achieve a better visu- 
alization of the system and often allow to avoid unessential technichal complications, even 
though one has to remain aware of potential oversimplifications. One good example is the 
problem of stability and order reduction (via BPS equations) of the nonlinear sigma and 
Skyrme models, which was considered in three dimensions by [3], leading to the discovery 
of instanton solutions. 

More specifically, there exists a 2+1 dimensional version of the Skyrme model which was 
originally introduced in 14], [l5| . and further investigated, e.g., in Jig], 17j], 18[ (for more 
recent results see, e.g., 19| - |22[). In this so-called baby Skyrme model, the dimensions 
both of the base space and of the target space are reduced by one, such that a topological 
degree characterizing field configurations and topological soliton solutions exist, again, in 
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close analogy to the full Skyrme model. This model has already found some independent 
physical applications, e.g., in condensed matter systems [23], or in brane cosmology where 
the solitons of the model induce co-dimension two branes [24]. The lagrangians of the Skyrme 
and baby Skyrme models are very similar, as well. In both models there exists a nonlinear 
sigma-model type kinetic term and a Skyrme term which is quartic in first derivatives. In 
addition, the potential term is mandatory in the baby Skyrme model for soliton solutions 
to exist, as a consequence of the Derrick scaling argument. The baby Skyrme model, too, 
has a natural coupling to the electromagnetic field. It turns out that the resulting baby 
Skyrme-Maxwell system in one lower dimension does allow for a nonlinear treatment with 
a reasonable numerical effort. Indeed, soliton solutions of the full coupled system have 
been calculated in [25] . and their electric and magnetic properties have been studied, too. 
Further, it turns out that in the gauged system the Skyrme term is not mandatory because 
the Maxwell term shows the same scaling behaviour, and the resulting "gauged nonlinear 
sigma model with a potential" obtained by skipping the Skyrme term has soliton solutions, 
too, see [26] . Both the gauged Skyrme model and the gauged nonlinear sigma model with 
potential have a BPS bound for the energy in terms of the topological degree, but nontrivial 
soliton solutions, in general, do not saturate this bound. Still, the gauged nonlinear sigma 
model does have genuine BPS soliton solutions for a specific choice of the potential 26] (see 
Section II). Apparently, the progress in the understanding of these 2+1 dimensional models 
is in part related to the reduced number of dimensions and has not yet provided us with 
all the necessary tools for a deeper analysis of the full 3+1 dimensional Skyrme-Maxwell 
system. So one might wonder whether there exists some further structure which is shared 
by the baby and the full Skyrme theory and which might bring us closer to this final goal. 
It is the purpose of the present article to give an affirmative answer to this question, but 
before presenting our proposal it is necessary to briefly review some recent results on BPS 
submodels of the Skyrme and baby Skyrme models. 

The field space (or target space) of the baby Skyrme model is given by the two-sphere 
S 2 . If we require that the Lagrangian density for this field variable depends on the fields 
and their first derivatives, is Poincare invariant and at most quadratic in time derivatives 
- in order to allow for a standard hamiltonian formulation - then the Lagrangian of the 
baby Skyrme model is the most general possibility. Specifically, the Skyrme term which is 
quartic in derivatives and is, at the same time, the square of the topological current density, 
already has the highest possible power of derivatives. The simplest submodel supporting 
topological solitons is the scale-invariant nonlinear sigma model [3] consisting only of the 
quadratic term, but otherwise both the potential and the quartic Skyrme term are required 
to maintain stability under Derrick scaling. The quadratic (sigma model) term, on the 
other hand, is not mandatory from this point of view, and the restricted model containing 
only the potential and the Skyrme term not only supports solitons but, additionally, has 
some further structure. It was first considered in 27| where the infinitely many base space 
symmetries of the static energy functional were found, and both exact static soliton solutions 
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and exact time-dependent topological Q-ball solutions were calculated. Then, in [28| it was 
found that this restricted baby Skyrme model has both a BPS bound and BPS soliton 
solutions saturating this bound. Further, it was demonstrated that the model has infinitely 
many target space Noether symmetries with their corresponding conservation laws (it has, 
in fact the zero curvature representation of generalized integrability 29j). A more geometric 
interpretation of the BPS bound and the corresponding BPS equations was given in (3oj . 

If the same conditions as above are imposed for the S 3 (or SU(2)) field space of the 
Skyrme model, then the most general lagrangian density consists of a potential term, the 
quadratic nonlinear sigma model term, the quartic Skyrme term and a sextic term which 
is, at the same time, the square of the topological current density. That is to say, the 
quartic Skyrme term of the baby Skyrme model has two possible generalizations which both 
are equally acceptable from the point of view of the Derrick scaling argument. The full 
generalized Skyrme model with all allowed terms present has been considered and applied 
to nucleons and nuclei (see e.g. 31], 33 ] ) , but a systematic search of higher solitons, 

analogous to the case of the standard Skyrme model, has not yet been performed to the 
best of our knowledge. This generalized Skyrme model, again, has a submodel (the "BPS 
Skyrme model" ) consisting of a potential and the sextic term, which shows all the additional 
structure mentioned above 34|, 35|]. The static energy functional has the volume-preserving 
diffeomorphism on base space as symmetries, which is interesting from the point of view of 
nuclear physics (nuclear liquid drop model), because these are precisely the symmetries 
of an incompressible liquid. Further, the BPS Skyrme model is integrable in the sense of 
generalized integrability and has an infinite number of target space Noether symmetries and 
the corresponding conservation laws. Thirdly, the BPS Skyrme model has a BPS bound 
and soliton solutions saturating this bound. This fact is of special importance because it 
resolves one of the main problems in the application of the Skyrme model to nuclear physics, 
namely the too large binding energies of nuclei which the standard Skyrme model predicts. 
As a consequence of the BPS property, classical soliton solutions in the BPS Skyrme model 
have zero binding energies, and realistic small values for the binding energies of nuclei may 
be achieved both by quantum corrections and by small contributions of additional terms in 
the lagrangian {36|, 56]. 

In view of the existence and of the rich structure of the integrable BPS Skyrme submodels 
described above, both in 2+1 and in 3+1 dimensions, one rather obvious question is whether 
these submodels may be gauged and whether the resulting gauged BPS Skyrme models 
maintain the integrability, the BPS property, and the close relationship between the 2+1 
dimensional and the 3+1 dimenional case. It is the purpose of the present article to study the 
2+1 dimensional case and to demonstrate that the resulting gauged BPS baby Skyrme model 
in 2+1 dimensions still is integrable and has both a BPS bound and BPS solitons which, in 
some cases, may even be calculated exactly. An investigation of the 3+1 dimensional case, 
following the lines developed in the present article, will be presented elsewhere. 

The paper is organized as follows. In Section II, we introduce the gauged baby Skyrme 
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model and briefly recapitulate some known results. Then we restrict to the gauged BPS 
baby Skyrme model which we want to investigate, derive its Euler-Lagrange equations 
and discuss the boundary value problem relevant for soliton solutions. In Section III we 
demonstrate that integrability still holds for the gauged model. In Section IV we derive 
the BPS property of the gauged BPS baby Skyrme model. Concretely, we derive the BPS 
bound in Section IV. A, whereas in Section IV.B we demonstrate that the BPS equations 
imply the static Euler-Lagrange equations. Section IV. C is dedicated to the issue of BPS 
soliton solutions. In Section V we present the numerical calculations of soliton solutions for 
several potentials and coupling constants. The soliton energies saturate the BPS bound in 
all cases. In Section VI we present exact BPS soliton solutions for some specific choices of 
the potential. Finally, Section VII contains our conclusions. We decided to give a rather 
detailed presentation of the model and its properties in this article, because some of the 
methods used and some of the results presented are quite different from known ones (e.g. 
the BPS bound is markedly different from all BPS bounds known to us), and a detailed 
understanding of the new features of this model and the new analytical methods will be 
indispensable both for the study of the analogous case in 3+1 dimensions and for its use in 
other contexts (e.g., to understand whether the new type of BPS bound presented here may 
be employed for other field theories). 



II. THE GAUGED MODEL 



The degrees of freedom of the ungauged baby Skyrme model are described by a three- 
component vector of scalar fields = (0i, 4>2,<t>z) with unit length, <fi 2 = <p\ + 0| + <f) 2 = 
1. That is to say, its target space manifold is a two-sphere S 2 . The Lagrangian of the 
(ungauged) standard baby Skyrme model consists of three terms, namely the nonlinear 
sigma model term L 2 quadratic in derivatives, the Skyrme term L 4 (which is quartic in 
first derivatives), and the potential term with lagrangian density C — — ^oV(<f>)- Here, the 
potential V((f>) always is a real, nonnegative function of its arguments. Both L 2 and L 4 are 
invariant under the full group of target space rotations SO (3) (i.e., rotations acting on 0), 
whereas the potential necessarily breaks part of this symmetry. Throughout this paper, we 
shall be concerned only with potentials which preserve a U(l) subgroup of this target space 
symmetry and which, further, have a unique vacuum (i.e., there is no spontaneous symmetry 
breaking). Concretely, we assume that 

V = V(n-$), V(1) = 0, (1) 

— * 

which implies that V is invariant w.r.t. rotations about the axis n and that = ft = const, 
is the (unique) vacuum configuration. In concrete calculations we shall always assume 



ft= (0,0,1) V = V(cf> 3 ) 



(2) 
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such that the vacuum value is 3 = (j)^ V ac = 1, but in more general expressions it will be useful 
to maintain the general vacuum vector n. Finite energy field configurations 4>(t, x) have to 
approach the vacuum configuration in the limit of large \x\ independent of the direction of 
x, lim^i^oo (f>(t, x) = n for all times. This allows for the one-point compactification of the 
physical base space R 2 which makes it topologically equivalent to the two-sphere S 2 . As a 
consequence, finite energy field configurations may be interpreted as maps S 2 — > S 2 which 
are characterized by an integer winding number (or topological degree) 

deg[0] = — / d 2 x4> ■ d x $ x d 2 <p = k , k e Z. (3) 
4vr J 

Due to its integer-valuedness this winding number is obviously conserved in time. 



A. The gauged baby Skyrme model 

The presence of the unbroken subgroup U(l) opens the possibility to gauge this subgroup 
by a U(l) gauge field whose dynamics is governed by the usual Maxwell term. The right 
coupling of the gauge field to the Skyrme field <fi is achieved by replacing ordinary partial 
derivatives with the covariant derivatives 

D a cj) = d a (f) + A a n x (4) 



as may be checked easily [26|, |25|. The resulting Lagrangian density (i.e., the gauged baby 
Skyrme model) is 

C = C 2 + U + £ + £ m (5) 
where C 2 is the gauged sigma model term 

C 2 = y(Di) 2 , (6) 



£4 is the Skyrme term 
and £ is the potential 



C, = -^(DjxD^) 2 , (7) 



£ o = -A o V(n-0). (8) 
Finally, Cm is the usual Maxwell term 

Cm = — ~^~~Fa/3 5 Fap — d a A/3 — dpA a . (9) 

Here, all fields are defined on 2+1 dimensional Minkowski space, and we use the metric with 
signature (+, — , — ). Further, the k = 0,2,4, M are non-negative dimensionful coupling 
constants. Following 25], we now introduce new coupling constants by extracting a common 
energy scale E from all terms in the Lagrangian. The resulting Lagrangian is 

L = E J d 2 x (^{Djf - ^(Dj x D^f - fi 2 V(n ■ $) - j^^) • (10) 
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Here, E has the dimension of energy and sets the energy scale of our model. In concrete 
calculations we shall always express energies in units of Eq, which is equivalent to setting 
Eq = 1. v is a dimensionless constant which takes the value v = 1 in the original gauged 
baby Skyrme model of |25| (for a nonzero u, v — 1 can always be achieved by an appropiate 
choice of the energy scale Eq), whereas we shall choose v — for the case of our restricted (or 
BPS) gauged baby Skyrme model. Further, A has the dimension of length, whereas /i and g 
have the dimensions of inverse length, g is, in fact, the electromagnetic coupling constant. 
In this paper we are mainly interested in static (soliton) solutions, so let us restrict to the 



static energy functional in what follows. Further, it has been demonstrated already in [25 1 
that static finite energy solutions must have zero electric field, at least for the spherically 
symmetric ansatz considered there. We shall assume this in our paper, i.e., we choose for 
the static gauge field 

A a (x) = (0,A 1 (x),A 2 {x)) , E 1 = E 2 = 0, B(x) = d x A 2 - d 2 A 1 (11) 

where B is the magnetic field in 2+1 dimensions and the electric field E^ is zero. Then the 
static energy functional which defines our static variational problem is 

E = E v (v, A, /i, g) = ^E j d 2 x (V(A0) 2 + X\Dj x Dj) 2 + 2^ 2 V{n ■ $) + ^B 2 ^j . 

(12) 



B. Some known energy bounds 

Before starting our calculations, it is useful to review some known results. It is well-known 
that the ungauged pure nonlinear sigma model with energy functional 

E u (is, A = 0, n = 0) = ^E a v 2 J d 2 x{dS) 2 (13) 

(here the subindex u stands for "ungauged", i.e., zero gauge field A a = 0) has a topological 
lower energy bound for static configurations 13] 



E u (u, A = 0,/i = 0) > 47r£ o z/ 2 |deg[0]|. (14) 

Further, the pure nonlinear sigma model has the meromorphic functions as static solutions 
saturating this bound (i.e., BPS solutions; in fact, all static finite energy solutions are BPS 
solutions). The same bound remains true for the full ungauged baby Skyrme model, 

^(z/,A, /U )>4vr J E; o z/ 2 |deg[0]| (15) 

because of the obvious inequality E^(u, A,/i) > E u {y, A = 0, \i — 0). Nontrivial solutions of 
the full baby Skyrme model do not saturate this bound. There exists an additional, more 
stringent bound for the full, ungauged baby Skyrme model which derives from the fact that 
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there exists a second, independent bound for the baby Skyrme model without the sigma 
model term, i.e., for the BPS baby Skyrme model E^(y = 0, A,//) (see Section IV). Again, 
nontrivial soliton solutions of the full baby Skyrme model in general do not saturate this 
bound. 

For the gauged baby Skyrme model, there exist some known BPS bounds, as well. First 
of all, for the gauged model without the quartic Skyrme term (i.e., for A = 0) there exists a 
specific choice of potential and coupling constants such that the resulting "gauged nonlinear 
sigma model with potential" has both a BPS bound and BPS solutions saturating this bound 



26|. The specific potential is 



K(n-0) = i(l-n-0) 2 (16) 

and the coupling constants must be chosen equal, \i = g (in addition to A = 0; further we 
set v — 1 without loss of generality). The resulting energy is 

E v '(y = \,\ = Q,n,g = n) = E v °(v = 1, A = 0, fi = l,g = 1) 

= ~E j d 2 x ((A0) 2 + 2V s (n • 0) + 5 2 ) (17) 

where we transformed \i and g = /i into /x = g = 1 by a variable transformation x — > fix to 
dimensionless new variables x. It was proved in [26j that this model obeys the same BPS 
bound like the ungauged pure nonlinear sigma model, 

E v '{v = 1, A = 0,fj, = l,g = 1) > AnE \deg[$]\ (18) 

and that there exist BPS solutions saturating the bound for winding number |deg[0]| > 1. 
For |deg[0]| = 1, on the other hand, finite energy solutions do not exist. 

The above BPS bound for the specific potential V s may be used, under certain condi- 
tions, to derive energy bounds for other potentials. To show it, we first rewrite the energy 
expression like 

E v (v = 1, A = 0, (i, g) = E J d 2 x ((A0) 2 + 2^V(n ■ $) + B 2 ^j (19) 

2 

where we transformed again to new dimensionless space coordinates x — > gx. If may 

2 

be bound by multiples of V s , i.e., > cyV s , then the above energy may be bound either 
by 47i\E o |deg[0] | (if cy > 1) or by 47ri?ocy|deg[0]| (if cy < 1). Let us choose the so-called 
"old baby Skyrme potential" 

V (n -$) = l-ft-$ (20) 
as an example. As a consequence of the inequality 
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we find the following energy bound 

^ > 1 : E v °{v = 1, A = 0,n,g) > An E \deg[$}\ 

2 

^<1: E v °(v = l,\ = Q,^g)>47rE ^\deg[$]\. (22) 
9 9 

In these models, however, soliton solutions will in general not saturate the bound, i.e., they 
are not of the BPS type. 

These bounds continue to hold for the full gauged baby Skyrme model, i.e., for A ^ 0, as 
a consequence of the obvious inequality E y (v = 1, A, /i, g) > E y [y — 1, A = 0, /i, g). For the 
old baby Skyrme model, e.g., we have (see 251 ]) 

^>1: E v °(v= l,X,fi,g) >47r£ o |deg[0]| 

2 

^<1: E v "(u = l,\,^g)>AnE ^\deg[$}\. (23) 
Again, these bounds are not saturated for nontrivial soliton solutions. 



C. The gauged BPS model 

Now we shall restrict to the case we finally want to discuss in detail, i.e., to the case 
without the sigma model term, v = 0, with lagrangian density 



c = _ T x ~ ^ v{n ■ ^ _ 4?^- (24) 



The Euler-Lagrange (EL) equations are derived by varying w.r.t. the fields and their deriva- 
tives, as usual. For the Skyrme field <fi there exists a minor subtlety related to the constraint 
2 = 1. This constraint should in principle be implemented by adding a Lagrange multiplyer 
term \(x)(l — <p 2 ) to the lagrangian density (we suppressed this term for the sake of brevity). 
As a consequence, the component of the EL equation in the direction of <ft is not dynamical 
but serves, instead, just to determine the Lagrange multiplier A(x). We may project on the 
dynamical part of the EL equation by multiplying the full EL equation with 0x , which leads 
to 

D Q K. a = -/j 2 n x 4>V'. (25) 
Further, the inhomogeneous Maxwell equation is 

d a F a ? = g 2 n-KP, (26) 

where 

t a = \ 2 D <P \$ ■ {D a 4> x D^) 1 . (27) 
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The corresponding energy functional is 



E = - / d 2 x (— (d $ x A0) + -E} + A 2 x D 2 pj + 2/i 2 \/ + —B' 

J \ 9 9 



(28) 



Next we assume n = (0, 0, 1) and the standard static ansatz 
0(r, 4>) 



( sin /(r) cosn0 
sin/(r) sinn^ I , A = A r = 0, = na(r) 



(29) 



\ cos/(r) 



then the electric field vanishes identically and B = na ^ . Further, this field configuration 
has winding number deg[0] = n provided that / obeys the appropriate boundary conditions, 
see below. The field equations can be reduced to 



V'(l + «) 2 sm 2 / + ^ 



2a' 



+ JL- sin fV' = 0(30) 



l + a\ 1 + a . 2 f . 

sin /H (1 + a) sin /cos/ 

r 



7?/ 



! A 2 



a"- 



AV(l + a)sin 2 // 



/2 



(31) 



where now V = V((f>s) = V(cos/) and V = V^ 3 . Further, we introduce the new variable 



to arrive at the following system of autonomous second order equations 



sin / \d y [f y {\ + a) 2 sin/] + J^'} = 
a yy = AV(1 + a) sin 2 // 2 . 



With the help of the function 



03 = cos/ = 1 - 2h a=-(l-cos/), /i y = - sin //„ 



this may be further simplified to 



sin f \ d y [h y (l + a) 2 ] 



4n 2 A 2 

2„2/i i „\^i2 



K > = 



(32) 

(33) 
(34) 

(35) 

(36) 
(37) 



<hm = A V(1 + a)4/ij 

where now V" = V(/i) and Vh = —IV. This system of two second order equations has, in 
general, four integration constants, i.e., a four dimensional solution manifold. We shall see 
in a moment that, on the other hand, a soliton solution imposes four boundary conditions, 
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therefore generically we expect at most one soliton solution for this spherically symmetric 
ansatz and for a given n and a given potential. At y = (i.e., r = 0), a soliton solution has 
to obey the two conditions 

h(0) = 1 /(0) = 7i , o(0) = 0. (38) 

The precise form of the further boundary conditions depends on the way the fields approach 
their vacuum values. If the soliton takes its vacuum value already at a finite radius r = r$ 
i.e., y = yn) and remains in the vacuum for y > y it is said to be of the compacton type 
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40| . 28]. In this case, the additional boundary conditions are 

a y (y = y ) = , h(y ) = h y (y ) = 0. (39) 

Apparently these are three more conditions, but we introduced the additional free constant 
yo (the compacton radius), so these correspond to two more conditions, and the total number 
of boundary conditions is four. Alternatively, if the field h approaches the vacuum expo- 
nentially (like exp —cy) or if the approach to the vacuum is power-like (like y~ c ) (here c is 
some positive real constant) then the fields reach their vacuum values in the limit y — > oo, 
so the boundary conditions for a soliton are 

lim %) = , lim a y (y) = (40) 

y— >oo y— ¥oo 

whereas \im y ^ 00 h y (y) = is not an independent condition but, rather, a consequence of 
the asymptotic form of the above field equations. Again, the total number of boundary 
conditions is four. 

We remark that in the ungauged model there is a simple relation between the potential 
V((p 3 ) and the approach to the vacuum, independent of the coupling constants. Indeed, if 
V is less than of a fourth power in small flucuations about the vacuum, then a possible 
soliton solution of the ungauged model must be of the compacton type, i.e., it takes its 
vacuum value already at a finite radius r = r (i.e., y = y ) and remains in the vacuum 
for y > y , [2Sjj]. Further, if the potential is exactly quartic in small fluctuations, then the 
approach of a soliton to its vacuum value is exponential, whereas if V goes to zero with a 
higher than fourth power then a soliton of the ungauged model approaches the vacuum in 
a power-like way. Observe that h itself is already quadratic in small fluctuations because 
h = (1/2)(1 - 3 ) = (1/2)(1 - y/l-(j>\- 0f) ~ (l/4)(0? + Therefore, a potential 
which is "quartic in small fluctuations" means a potential which is "quadratic in h" for 
small h. The "old" baby Skyrme potential (1201) . e.g., is quadratic in small fluctuations, 
whereas the potential (TIBjl is quartic. At this point the obvious question shows up whether 
the simple relation between the potential and the approach to the vacuum continues to 
hold in the gauged model, or whether the approach to the vacuum will depend both on the 
potential and on the values of the coupling constants. For the specific case of the so-called 
old baby Skyrme potential it is easy to see that the potential alone is sufficient to determine 
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the asymptotics, see next subsection. For the general case it is more difficult to answer 
this question, but we will find strong theoretical indications that still the potential alone is 
sufficient to determine the asymptotic behaviour, see Section VI. 



D. The old baby Skyrme potential 
Although we will consider different potentials, most of our numerical calculations will be 



done for the old baby Skyrme potential (1201) . One reason is that this is the case considered 
in 25J]. A further reason is related to the fact that for the old baby Skyrme potential the 
system of second order equations (1361) . (|37|) has a simple first integral. Indeed, for the old 
baby Skyrme potential these two equations read 

sill /{^Nl+«f]-^}=0 (41) 
a yy = AV(1 + a)Ah 2 y . (42) 

The first equation is solved by 

sin/ = f = h = (43) 
or by setting the expression in brackets equal to zero, which can be easily integrated to give 



h y (l + a) 2 = ^- 2 (y-y ). (44) 

It already follows from this first integral that y^ is the compacton radius and that the 
nontrivial solution of the above equation for y < y must be joined with the vacuum solution 
h = for y > yo if we want to avoid a solution which grows indefinitely for large y. Formally, 
this first integral may be further integrated to give 

The magnetic function a can be obtained from eq. (1421) . where the baby skyrmion profile 
function is expressed in terms of a by (JHJ). Then, we get 

(l + a)\ y = X 2 g 2 (-^\ (y-y ) 2 (46) 

or 

(1 + aYa yy = (47) 
outside of the compact baby skyrmion. It gives a trivial solution outside the compacton 

a{r > r ) = clr = const. (48) 
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This implies that also the magnetic field is confined inside the compact baby skyrmion. 
Hence, we still have genuinely compact objects which interact only via a contact interactions. 
Obviously we may also construct multisoliton configurations exactly as in the standard BPS 
baby Skyrme model. The equation for a can be simplified by introducing b(y) = 1 + a(y) 

and = AV (^) 2 - Further z = f3 l '\y - y Q ), then 

b% z = z\ (49) 

Unfortunately, we were not able to find analytic solutions to this equation, so some numerics 
is still required. Of course, a solution of the magnetic equation should be inserted into the 
first order equation ( ]44"1) . where we also have to implement the boundary conditions for a 
compacton solution. Numerical solutions will be calculated and described in detail in section 
V. 

We remark that the possibility to partially integrate the static field equations already 
points towards the possible existence of a BPS bound and BPS equations for the static 
system. We shall see later, in section IV, that this is indeed the case, i.e., there exists a 
BPS bound, and soliton solutions are, in fact, solutions to the related BPS equations. The 
existence of such BPS solutions will, however, be related to some nontrivial conditions for 
the potential which not all potentials satisfy. As a consequence, there exist potentials such 
that the ungauged models has BPS soliton solutions, whereas the gauged model with the 
same potential does not support soliton solutions, even for arbitrarily small electromagnetic 
coupling constant g. 



III. INTEGRABILITY, SYMMETRIES AND CONSERVATION LAWS 



The ungauged model 



C 



<f>xdp<f>) -fi 2 V(n 



(50) 



has an infinite number of target space symmetries with their corresponding Noether currents 
and conservation laws The model possesses, in fact, the zero curvature representation 
of generalized integrability j^]. It is not difficult to understand the geometric origin of 
these symmetries. The quartic Skyrme term alone is invariant under the full group of 
area-preserving diffeomorphisms of the target space S 2 , because it is just the square of the 
pullback of the corresponding area two-form 28|]. The potential is only invariant under the 
subgroup which does not change n ■ (ft = 3 , i.e., under the abelian subgroup [41| 



ii -> 4>'x = cos/(0 3 )0i -sin/(0 3 )0 2 
k -» 02 = sin/(0 3 )0i + cos/(0 3 )0 2 

!>3 0' = 03 



(51) 
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where f{4>$) is an arbitrary function of its argument, which makes the symmetry group still 
infinite dimensional. It turns out that the gauged model maintains exactly the same abelian 
symmetry group and the corresponding conserved currents. 

Remark: It might appear that in the gauged model the above transformations are just a 
subset of the full group of gauge transfomations and, therefore, should not provide conser- 
vation laws, because it is a well-known fact that gauge transformations do not give rise to 
nontrivial conservation laws. This is, however, not true. The important point is that the 
gauged model is separately invariant under the above transformations (1511 and under the 
transformations 

A a -> A' a = A a - d a g{<f> 3 ) (52) 

whereas only the joint transformations fl5"Tj) and fl52l) with / = g are gauge transformations. 
It follows that the transformations (ISTjl are genuine symmetry transformations which give 
rise to conserved Noether currents and conserved charges. As a consequence of the above 
argument, the symmetry transformations (1521) provide exactly (minus) the same conserved 
charges. 

For an explicit calculation of the conserved currents using the methods of generalized 
integrability it is useful to switch to the CP 1 formulation of the gauged BPS baby Skyrme 
model. The Lagrangian density reads 

C = A 2 1 (D a uD(su* - D a u*D p uf - fi 2 V(uu*) + -^i^ (53) 
(1 + \u\ J 4(7 

where the complex scalar field u is related to the unit vector via stereographic projection, 

4> = - — ; — — (u + u, —i(u — u), \u\ 2 — l) (54) 

and the covariant derivatives are 

D a u = u a — ieA a u, D a u* = u* a + ieA a u*. (55) 
The currents resulting from generalized integrability read 

j a = iG'(u*Ti* a - UTC a ) (56) 

where G = G(uu*) is an arbitrary function of its argument and ir a is the canonical momen- 
tum 

- D p u*D a u)D^u (57) 

- D p u*D a u)D^u*. (58) 

u*Dpuf + /2 2 V'u\ (59) 



, 8C 4A 2 , 

du* a (1 + | M | 2 ) 4V p 



dC 4A 2 

du« ~ ~(1 + \u\ 2 ) 

The field equations are 



4u* 

D ^ a = X \, , i [D a uD p u* - D c 
(1 + \urr 
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An i 

D » n * a = A ' n , i 12^5 ( D * uD P u * ~ D a u*D p u) 2 + n 2 V'u. (60) 

Then, 

,9 a j a = iG'(u* a 7T* a - M Q 7T Q + U*O a 7l* a - m9 q 7T Q ) + iG"(uU* + U*U a )(u*7T* a - U7T a ). (61) 

The first parenthesis is 

... = <7r* a - u a 7c a + ieA a u*7r* a + ieA aU 7c a - fi 2 V'uu* + fi 2 V'uu*+ (62) 

+ A 2 (1 ^* 2)5 (D a uDpu* - D a u*Dpuf - A 2 (1 ^* 2)5 (D a uDpu* - D a u*Dpu) 2 (63) 

4A 2 

= D a u*n* a -D a un a = —^-(DpuD^-D^D^D^uD^+D^D^u*) = (64) 

(1 + \u\ ) 

as we contract an antisymmetric tensor with a symmetric one. Further, the second term 
leads to 

4A 2 

... = iG"(uu* a + u*u a )- —-{D p uD a u* - D l3 u*D a u)( y u*D /3 u + uDpu*) (65) 

(1 + \u\ ) 

4A 2 

= iG "i r^Kj( Dl3uDau * ~ D l3 u f D a u){uu* a + u*u a ){uu* B + u*ur) = 0. (66) 

(1 + \u\ 2 ) 4 p 

Therefore, the current is conserved for arbitrary G(uu*), and there exists an infinite number 

of conservation laws, as announced. 

The static energy functional of the ungauged model has the group of area-preserving 

diffeomorphisms on base space as an additional group of (non-Noether) symmetries. Let us 

now demonstrate that this symmetry, too, is maintained in the gauged model. First of all, the 

Skyrme term of the gauged model is obtained from the Skyrme term of the ungauged model 

by simply replacing the partial derivatives dj by covariant derivatives Dj. Further, dj and 

Dj have the same behaviour under coordinate transformations (both are covectors), so the 

transformation of the Skyrme term under general coordinate transformations (and, therefore, 

also under the subgroup of are-preserving diffeomorphisms) is the same for the gauged and 

the ungauged case. For the potential, obviously nothing changes. The potential itself is a 

scalar and is, therefore, invariant under general coordinate transformations. The invariance 

of the static energy functional under area-preserving diffeomorphisms follows, therefore, 

from the same invariance of the area two-form d 2 x. Finally, the gauged model contains the 

Maxwell term. But for a static, purely magnetic configuration = (0, Ai(x), A 2 (x)), the 

Maxwell term is proportional to B 2 where B is the magnetic field and transforms like a 

pseudoscalar under coordinate transformations on the two-dimensional base space of static 

configurations. It is, therefore, invariant under orientation-preserving diffeomorphisms and, 

consequently, under area-preserving diffeomorphisms. We conclude that the static energy 

functional of the gauged model, like the ungauged one, is invariant under area-preserving 

base space diffeomorphisms and shows the same degeneracy among field configurations with 

different shapes but the same area. 
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IV. THE BPS BOUND 



The BPS bound of the ungauged model was first found in [42j (as a contribution to an 
improved bound for the full baby Skyrme model), whereas the BPS nature of the restricted 
model was proved in [3] and in 30]. Here we follow the geometric discussion of 30]. In a 
first step, we want to briefly recall the BPS bound of the ungauged model, because we will 
need this result later on. The static energy functional of the ungauged model is 



E 



-En 

2 

— En 

2 



d x 



X 2 ldKf> x d 2 <pj + 2fx 2 V((j) ■ n) 
d 2 x [AY + 2/iV(0 ; 



(67) 
(68) 



where 



q = 4> ■ di<p x d 2 <f) (69) 
(i.e., q is An times the topological charge density). This energy functional leads to the bound 

E = ^E J d 2 x (\q ± fiVWy T E \fi J d 2 xq\[W > TE \fi J d 2 xq\fW (70) 

with equality if the BPS equation 

Xq±fiV2V = (71) 

is satisfied. It remains to show that the bound is, in fact, topological (i.e. does not depend 
on the field configuration <fi). This follows from the following fact. As said already, the 
Skyrme field <ft defines a map $ from the one-point compactified base space M. 2 to the target 
space S 2 , 3> : 1R 2 — > S 2 . Further, the two-form d 2 xq is just the pullback under this map of 
the area two-form Q on the target space S 2 , i.e., 



d 2 xq = 



(72) 



It follows that J d 2 xq\/2V is just An (i.e., the area of the two-sphere) times the average 
value of \/2V on target space times the times covers the target space while x covers the 
base space once (i.e., the topological degree or winding number), that is 



d 2 xq^2V( ( j )3 ) = Andeg[$]{VW) S 2. 



(73) 



Obviously, this expression does not depend on the field configuration (j)(x). For a given 
potential ^(^3), the average value (y/2V)s^ takes a fixed, given value, so this expression only 
depends on the winding number deg [</>], i.e., it is a topological quantity. The corresponding 
energy bound reads 

E >An\fxE \deg[(p}(V2V) S 2\. (74) 

Further, there exist BPS soliton solutions which satisfy the BPS equation ( TfTj) and, therefore, 
saturate this BPS bound, see [28|, [30 ] . 
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A. The BPS bound 



The static energy functional of the gauged model is 



E = ^E I <f v 



X 2 [D l( px D 2 <P) +2^V(<j)-n) + —B' 



1 



E / d 2 x 



X 2 Q 2 + 2pt 2 V (0 3 ) + - 2 B< 

g 2 



where 



Q = <f> ■ D 1( f) x D 2( f). 

It further holds that 

— # 

Q = q + eijAidjin ■ 0). 

From now on we choose n — (0, 0, 1), n ■ $ — <p 3 . 

Next we consider the following non-negative expression 



(75) 
(76) 

(77) 
(78) 



^E I d 2 x 



X 2 {Q-w{4> 3 )f + -AB + b{<t, 3 )f 



r 



(79) 



where w(0 3 ) and b(<f) 3 ) are functions of 3 which we shall determine below. The non- negative 
expression may be written like 



^En / d 2 x 



\ 2 Q 2 + \ 2 w 2 + j 2 b 2 + j 2 B 2 - 2\ 2 wq - 2X 2 we ij A i d j( p 3 + 2^e ij {d i A j )b 



1 



(80) 



where we used B = e^diAj. The last two terms combine into a total derivative if we choose 



r<t>3 

b( ( p 3 )=g 2 \ 2 W( ( p 3 )=g 2 \ 2 / dtwit). 

J <l>3,vac 



Indeed, the last two terms now give 



E J d 2 xX 2 e ij d j (A i W) 



(81) 



(82) 



and this gives zero because by construction W((fi 3 = (p 3jVac = 1) = 0. The non-negative 
expression therefore reads 



-£o / d 2 x 



X 2 Q 2 + X 2 W' 2 + g 2 X 4 W 2 + ^B 2 - 2X 2 W'q 



1 

? 



If we now require that W obeys 

X 2 W' 2 + g 2 X 4 W 2 = 2/j 2 V 
then we get the BPS bound for the energy 

E > E X 2 J d 2 xqW = 47r\k\E X 2 (W') S 2 



(83) 



(84) 



(85) 
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where k = deg[0] is the winding number. There is a sign ambiguity in the choice of W, so 
we may always choose the sign such that (W')sn > 0. 

We remark that eq. for W(4>3) is analogous to the "superpotential equation" which 
determines a superpotential W(<j)) in terms of apotential V(4>) for a real scalar field 6 in the 



context of selfgravitating domain walls 43[ - 48J or scalar field inflation models 49] - 52 ]. 



There the superpotential allows to reduce the domain wall or cosmological equations to a 
first order form. Further, the relation between potential and superpotential is completely 
equivalent to the relation in the corresponding (dimensionally reduced) supergravity theories, 
and the method is therefore called "fake supergravity". We shall find that in our case the 
auxiliary function W allows a reduction to a first order system, as well. The main difference 
is that in fake supergravity the two terms proportional to W 2 and W' 2 enter with different 
signs in the superpotential equation, whereas they enter with the same sign in our eq. ()84l) , 
which implies some qualitative differences in the solution space, as we shall see in a moment 
(we remark that in the context of extremal, supersymmetric black holes the superpotential 
equation shows up, as well, and there both terms enter with the same sign, as in our case, 
see e.g. 0|). In view of these similarities, we shall call W the "superpotential" and eq. 
fl84|) the "superpotential equation" in what follows. 

The formal energy bound (185j) only gives a genuine BPS bound provided that either W 
itself or at least (W')s^ are uniquely determined. On the other hand, eq. (184j) is a first 
order ODE, so it seems that it provides a one-parameter family of solutions W. It turns 
out, however, that this is not true, and eq. (|84|) only has one unique solution. The reason 
is as follows. If we want to find a local solution of eq. (184"]) which is valid in the vicinity 
of a point 3 = c, where —1 < c < 1, then, indeed, there exists a one-parameter family of 
solutions. We may choose the "initial value" W(c) from the interval 

Vmc) < W(c) < J^- Jw{c), (86) 



gX 2 w " ~~ g^ 2 

and each choice produces one solution. The important point is that we require a solution 
which exists in the whole interval — 1 < 3 < 1, specifically at the vacuum value 3 = 1 
where the potential is zero, V(l) = 0. But for this vacuum value, the "interval" of eq. (18"B"|) 
collapses to the point W(l) = 0, and this is the only allowed "initial value" at 3 = 1 for 
the differential equation (!H4|) . The solution is, therefore, uniquely determined in the vicinity 
of the vacuum 03 = 1. A different (and very important) question is whether this unique 
solution can be extended to the whole interval 3 G [—1, 1], i.e., to the whole target space, 
see Section IV.C Assuming this for the moment, we find for (W')ga 

p />2n pit 

An{W') S 2 = / dnW = / dip d0sin0W(cos0) 
Js 2 Jo Jo 

= 2tt J dtW'{t) = 2tt(W(1) - W(-l)) = -2ttW(-1) (87) 
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and for the BPS bound 

E > E X 2 J d 2 xqW = 27i\k\E X 2 \W(-l)\. (88) 

B. The BPS equations 

It still remains to prove that the two BPS equations imply the static field equations. The 
BPS equations are 

Q = W (89) 

B = -g 2 X 2 W, (90) 

where we used the condition that the BPS equations hold if and only if the non-negative 
expression (179"]) is zero, together with the expression ( IHTj) for 6(^3) and w = W. On the 
other hand, the static second order field equations are 

X 2 e lj D i [{D j $)Q] = -|iV'n x (91) 
3 t F lj = g 2 X 2 n- D k $((f- D j (f x D k $). (92) 



First we prove that the two BPS equations (1891) and (1901) imply the inhomogeneous static 
Maxwell equation (192]) . With n ■ D k <p = dk4>3 the Maxwell equation reads more explicitly 

d k B = -g 2 X 2 d k ^Q. (93) 

On the other hand, the partial derivative of the second BPS equation is 

8 k B = -g 2 X 2 W'd k <p 3 = -g 2 X 2 Qd k ^ (94) 

where we used the first BPS equation in the last step. It is, therefore, identical to the 
Maxwell equation. 

Finally we prove that the two BPS equations imply the static field equation (1911) for the 
Skyrme field. First we observe that Eq. ( !9T|) may be re-expressed like 

D 2 $d l Q - D x 4>d 2 Q + n x $BQ = -X~ 2 ii 2 V'n x $. (95) 

Then we use the equations 

fi 2 V' = X 2 W'W" + g 2 X i WW\ d k Q = W"d k {n ■ $) (96) 

and the second BPS equation to arrive at 

(D 2 $d 1 (n ■ 0) - D 1( fd 2 (n ■ $))W" - g 2 X 2 WW'n x = -(WW + g 2 X 2 WW')n x (97) 

or at 

D 2 <pd!(n ■ 0) - D 1( pd 2 (n ■ 0) = -W'ft x 0. (98) 
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Inserting for the covariant derivatives, this reads 

<9 2 0<9i(n • 0) - d 1( fd 2 (n ■ $) = n x 0[A 1( 9 2 (n ■ 0) - A 2 d 1 (n ■ 0) - W] (99) 

and using a last time the BPS equation W = Q and the expression f!78|) for Q we find that 
all terms containing the gauge field disappear and we are left with 

d 2 <f)di(n ■ 0) — di4>d 2 {n ■ 0) = — n x 0(0 • <9i0 x <9 2 0). (100) 

We still have to demonstrate tha this equation is, in fact, an identity. First, we observe that 
both sides in this equation are perpendicular to and to n, therefore we may project the 
only nontrivial equation by the scalar product with the vector n x which results in 

(n • x d 2 (pdi(n ■ 0) — (n • x di<f)d 2 (n ■ 0) = —(1 — (n ■ 0) 2 )0 • <9i0 x <9 2 0. (101) 

Finally, we need the identity 

n ■ <9i0 x <9 2 = {n ■ 0)0 • <9i0 x <9 2 (102) 

which follows from the fact that <9i0 and d 2 <p span a plane perpendicular to to arrive at 

(n ■ x d 2 (pdi(n ■ <p) — (n ■ <p x di<j)d 2 (n • 0) + ■ <9i0 x <9 2 — (n • 0)n ■ <9i0 x 9 2 = 0. (103) 

But this last equation is an identity as an immediate consequence of the Schouten identity 
in three dimensions, 

^ab^cde ^ac^deb _|_ ^ad^ebc ^ae^bcd g (104) 

Indeed, contracting the Schouten identity with n a n b (j) c di4> d d 2 (j) e we arrive at the above equa- 
tion. Therefore, both static field equations are consequences of the BPS equations, which is 
what we wanted to prove. 

C. BPS soliton solutions 

1. General remarks 

Eq. fl84|) has one unique solution defined by the initial condition 1^(1) = 0, W'(l) = 
at the vacuum value 03 = 1. If this initial condition can be integrated such that the 
solution covers the whole interval — 1 < 03 < 1, then there exists a unique, well-defined 
BPS bound for this model (i.e., for the corresponding potential with vacuum at 03 = 1 and 
for the corresponding choice of the coupling constants A, fi and g). It happens that for 
some potentials (or coupling constants) Eq. (1H4|) cannot be integrated for the whole interval 
— 1 < 03 < 1- In these cases, the BPS bound derived in this section does not apply for 
that model. We shall find both possibilities in the following. Concretely, both for the "old" 
baby Skyrme potential V Q = 1 — 3 = 2h and for the potential V D 2 = 2V S = 4h 2 we will 
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find that a solution to eq. (1841) exists in the whole interval 3 G [—1, 1], i.e., h G [0, 1]. For 
these potentials we find that soliton solutions exist and are, in fact, solutions to the BPS 
equations ( 1891) and ( 1901) . that is, they saturate the BPS bound ( IHHi) . A global solution to 
eq. seems to exist for the class of potentials V a = (2h) a (we found them by numerical 
integration for different values of a), and we believe that BPS soliton solutions exist for all 
values of a for which eq. has a global solution. For the so-called "new" baby Skyrme 
potential 

V n = ±(l-<f>l) = h(l-h), (105) 

on the other hand, eq. (184j) does not have a global solution. Starting at h — with 
the initial condition W{h = 0) = (and, of course, Wh(h = 0) = 0), the numerical 
integration develops a singularity somewhere between h = 1/2 and h = 1 (the precise 
position of the singularity depends on the values of the coupling constants A, /i and g). 
Further, we are not able to find numerical soliton solutions for the new baby Skyrme 
potential, neither BPS nor non-BPS ones. These numerical findings, and the analyti- 
cal results described below (see Sections IV.C.2, IV.C.3), motivate the following conjectures. 

Conjecture 1: If a given restricted gauged baby Skyrme model (i.e., with a given 
potential and fixed values of the coupling constants A, /i and g) has topological soliton 
solutions at all, then these solutions are BPS solutions, i.e., solutions of the two BPS 
equations (189|) . (f90|) . where the superpotential W is the solution of the superpotential 
equation with boundary condition W((ps = 1) = W(h = 0) = 0. 

Corollary 1: The global existence of the superpotential (i.e., the existence of a solution to 
Eq. in the whole interval 03 G [—1,1], i.e., h G [0, 1]) is a necessary condition for the 
existence of soliton solutions. 

Conjecture 2: The superpotential equation and its solution W completely determine 
the existence of soliton solutions. More specifically, in the generic case solitons exist if and 
only if the superpotential equation (IMj) has a global solution on the whole interval h G [0, 1] 
(i.e., on the whole target space) and this solution obeys W ^ for — 1 < 3 < 1 (i.e., 
W h ^ for < h < 1). 

Remark: by "generic" we mean that there might exist exceptions, i.e, potentials where 
the corresponding superpotential exists globally and obeys Wh = at some points but, 
nevertheless, supports BPS solitons. If they exist at all, these potentials will, however, be 
rare in the sense that they require some fine-tuning of parameters, see the discussion in the 
next subsection. 

Conjecture 3: Potentials which have one vacuum at 03 = 1 and are strictly monotonous 
in the open interval 03 G (—1, 1) (i.e., potentials which have one vacuum at h = and are 
strictly monotonous in the open interval h G (0, 1)) satisfy the conditions of Conjecture 2 
and support, therefore, topological BPS solitons. 

Remark: we do not think that the conditions on potentials in Conjecture 3 are necessary. 
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That is to say, most likely there exist potentials which do not satisfy these conditions and 
still support BPS solitons. 

A rigorous proof of these conjectures is probably quite difficult and is certainly beyond the 
scope of the present article. 



2. The problem with Wh = 



We demanded in conjecture 2 not only the global existence of the superpotential W but 
also the absence of local extrema in the open interval < h < 1, so let us explain the 
problems related to Wh = 0. We restrict our discussion to the spherically symmetric ansatz 
for the BPS equations, i.e., we assume that if soliton solutions exist, at all, then they should 
also exist in the spherically symmetric subsector. We use the function h = (1/2) (1 — 03 ) 
instead of 03, then the superpotential equation for W is 

^-W 2 + g 2 \ A W 2 = 2fi 2 V(h) (106) 

and the two BPS equations for the spherically symmetric ansatz read 

2nhy(l + a) = --W h (107) 

na y = -g 2 \ 2 W (108) 

where a(y) is the angular part of the gauge field and y = r 2 /2. Next we resolve eq. fll06j) 
for W h , 

W h = \I%^V -Ag 2 \ 2 W 2 (109) 



and calculate the derivative 



,2 



w hh = 4r (tj v >i ~ 9 2 ^ 2 ww h ) . (i in ) 



w h \x 

For Wh = this is nonsingular only provided that Vh = 0. If Whh is singular at a point h = h s 
where < h s < 1, then the integration breaks down at this point and cannot be extended 
further. The superpotential, therefore, does not exist globally, and the corresponding theory 
does not support BPS solitons. Further, this is the generic case in the sense that if the 
integration of the superpotential equation produces Wh = for some value h s , then even 
for potentials which obey Vh = for some values hi, generically h s will not coincide with 
any of the hi, and some fine-tuning of the parameters is required to make them coincide. 
Specifically, if Vh ^ in the whole open interval, then Wh = automatically produces a 
singularity. We remark, however, that numerically we found that for such potentials Wh = 
never occurs, which motivated our Conjecture 3. 

If Wh = occurs exactly at a point h s where Vh = 0, too, then the superpotential may 
exist globally (we shall see an example in Section VI) but, still, this does not imply that 
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BPS solitons exist. If Wh, e.g. has exactly one zero coinciding with one zero of Vh, then it 
follows from the first BPS equation (11071) that h y changes sign at the point where Wh = 
and is, therefore, positive either near h = or near h = 1. But this is incompatible with 
the restriction h 6 [0,1] together with the boundary conditions h(0) = 1, h(oo) = 0. It 
might appear that we could avoid this conclusion by assuming that it is the function 1 + a 
which changes sign at Wh = 0, but we shall see that this is impossible because a satisfies 
the inequality a(y) > — 1 V y, see Subsection IV. C. 6. It follows that all W where Wh has an 
odd number of zeros in the open interval < h < 1 are forbidden. 

This still leaves the possibility of globally existing superpotentials with an even number 
of zeros of Wh which support soliton solutions. At the moment we cannot exclude this 
possibility, but if it exists then it requires a large amount of fine-tuning. The derivative of 
the potential, Vh, must have at least the same number of zeros, and the parameters of the 
potential must be fine-tuned such that the positions of all the zeros of Wh coincide with the 
positions of (some or all of) the zeros of Vh- This is the fine-tuning mentioned in the remark 
after Conjecture 2. 

3. The boundary value problem of BPS solitons 

Each solution of the BPS equations (and of the superpotential equation (11061) ) is a solution 
of the static field equations. The converse, however, is not true in general. This is especially 
easy to see for the radially symmetric ansatz of ( l29l) . There, the reason is that each BPS 
equation provides one integration constant and eq. (I106[) provides none, therefore there is 
a total of two integration constants, and the space of solutions is two-dimensional. The 
original static field equations, on the other hand, provide four integration constants for the 
radially symmetric ansatz, so their solution space is four-dimensional. 
The important question is, of course, whether a soliton solution can be a solution of the BPS 
equations. A simple count of the boundary conditions which a topological soliton has to obey 
in the radially symmetric case seems to indicate that this is impossible, because a soliton 
solution has to fulfill the four boundary conditions fl38l) . fl39|) in the case of compactons (or 
fl38|) . (HI)]) for non-compact solitons), which requires four integration constants. 
What may still happen is that the condition W(h = 0) = 0, which uniquely fixes the solution 
of eq. (11061) . implies, at the same time, that the two boundary conditions at the compacton 
boundary (for compact solitons) or in the limit y — > oo (for non-compact solitons) are 
automatically satisfied, so that only two more boundary conditions at the center y = 
are left. In this case, the true soliton solution could be a BPS solution in the general 
case. Now we shall see that this is exactly what happens. We will explicitly discuss the 
case of a compact soliton, but the non-compact case is completely equivalent. Further, we 
shall restrict our discussion to strictly monotonous potentials where V(h = 0), Vh > for 
< h < 1 which requires Wh ^ in the same interval. 



24 



It obviously follows from the first BPS equation (I107P that h y = =>• Wh = 0. Further, 
we know that the unique solution of eq. (II 06ft obeys W(h = 0) = Wh(h = 0) = at the 
vacuum value h = 0. If h = is the only point where = is possible, then we can 
conclude that h y = =>■ ft = 0, which is exactly the first compacton boundary condition. 
But we know already that Wh = is forbidden in the open interval < h < 1, therefore 
W h = implies either h = 0,W = oi h = 1,W = (fi/g\ 2 )^2V(l). 

Further, we know that the condition W(0) = 0, Wh{0) = at h = always holds for the 
solution W of the superpotential equation fjl06j) . because it is our boundary condition for 
this unique solution. The condition Wh(l) = 0, W(l) = (/// ' g\ 2 )^2V{l) at h = 1, on 
the other hand, constitutes an additional boundary condition for this unique solution and 
will, therefore, hold only for exceptional, finetuned potentials. Specifically, it does not hold 
for the concrete potentials which we consider in this paper as may be checked easily by a 
numerical integration, therefore for these potentials we may conclude that h y = =>- h = 
which is exactly the first compacton boundary condition. 

Equation fl 106[) together with the two BPS equations (I107P and (11081) . therefore, automat- 
ically imply the first compacton boundary condition (3y such that h y (y ) = h(yo) = 0) 
in these cases. Further, h(y ) = =>- ^(^(^o)) = and the second BPS equation (I108P 
imply the second compacton boundary condition a y (yo) = 0. Therefore, the two integration 
constants of the two BPS equations will, in general, be sufficient to fulfill the two remain- 
ing boundary conditions h(0) = 1 and a(0) = 0. The compacton solutions are, therefore, 
BPS solutions, at least for the specific class of potentials considered in this subsection. Our 
numerical calculations completely confirm this result. 



We mentioned already that the specific potentials V a ~ h a all seem to allow for a global 
solution to the superpotential equation and, therefore, for a BPS bound. We show the result 
of a numerical integration for the values a = 1, i.e., the old baby Skyrme potential, and for 
a = 2, i.e., for the potential V s = 2h 2 , in Figures [1] and [2], respectively. If a BPS bound 
exists, then in the limit of vanishing electromagnetic coupling g the resulting BPS bound is 
just the BPS bound of the ungauged BPS baby Skyrme model, and there exists an exact 
expansion for Eq. ( 11061) in g 2 . We may, therefore, use this power series expansion to get 
exact (instead of just numerical) values for the BPS bound. Concretely, the power series 
expansion for W may be found from eq. (11091) . re-expressed as 



by the following three-step process, i) solve the equation iteratively for Wh, ii) expand the 
second root in a power series in g 2 and iii) integrate the resulting power series w.r.t. h, 



4- Small g expansion 




(111) 
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respecting the boundary condition W(h = 0) = 0. Let us consider the class of potentials 
V = h 2a , a > 0, as an example. Applying the procedure, we easily find that up to first 
order in g 2 the resulting W reads 

W W(h) = -p-^-h a+1 ( 1 - ^ -g 2 h 2 ) . (112) 

1 ; A(a + 1) V (a + l)(a + 3) y J 1 ; 

Evaluating this expression at h = 1 gives the BPS bound. Here, the leading order g° is 
the bound of the ungauged model, and the leading correction for small g is of order g 2 and 
negative. All these results, including the precise numerical values, are confirmed by our 
numerical calculations. On the other hand, there does not seem to exist an expansion of 
the superpotential equation for large g. Numerically, we find that the BPS bound (and, 
therefore, also the soliton energy) for large g behaves like g~ l in the case of the old baby 
Skyrme potential, see Section V. 

w, w 




FIG. 1: Solution of the superpotential equation (11061) for the old baby Skyrme potential 
V = 2h, for the coupling constants A = 1, g — 1 and /i 2 = 1/2 (W . . . dashed line). The 

solution exists in the whole interval h G [0, 1]. 



5. Potentials with two vacua 

Here we want to show that for potentials with two vacua BPS soliton solutions cannot 
exist. Concretely we shall focus on potentials with their two vacua at h = 0,1, although 
the generalization to other cases poses no difficulty. One specific example is given by the 
new baby Skyrme potential V n = h(l — h), but the same arguments apply for the general 
case. It is easy to understand why generically the superpotential equation does not have a 
global solution. The problem is that, as the potential vanishes at the two points h — and 
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FIG. 2: Solution of the superpotential equation HI 06[) for the potential V s = 2h 2 , for the 
coupling constants \ = 1, g = 1 and /z 2 = 1/2 (W . . . dashed line). The solution exists in 

the whole interval h G [0, 1]. 



h — 1, the putative solution would have to obey the two boundary conditions W(0) = 
and W(l) = 0, but these are too many conditions for a first order equation. We show the 
result of a numerical integration for the new baby Skyrme potential in Figure [31 where the 
singularity is clearly visible. As already mentioned in Subsection IV.C.l, for the new baby 
Skyrme potential there do not seem to exist soliton solutions, at all. It is not surprising that 
there are no BPS soliton solutions, because we need the superpotential in the BPS equations, 
so if there is no global superpotential solution we do not expect true (globally existing) BPS 
soliton solutions. On the other hand, it is not so obvious why there are no soliton solutions, 
at all (i.e., neither BPS nor non-BPS soliton solutions), although we found some indications 
that all soliton solutions are, in fact, BPS solutions (see the results of Subsection IV.C.3). 
The non-existence of solitons is in some sense surprising because the non-gauged BPS baby 
Skyrme model (i.e., the case g = 0) with the new baby Skyrme potential has a perfectly 
well-defined BPS bound and soliton solutions saturating this bound. Nevertheless, both the 
BPS bound and the soliton solutions cease to exist for arbitrarily small but nonzero gauge 
coupling g. The resolution of this puzzle resides, of course, in the superpotential equation 
which tells us that the superpotential W does not have to satisfy any boundary condition 
in the non-gauge case with g strictly zero. On the other hand, for arbitrarily small but 
nonzero g, W has to obey the two boundary conditions W(h = 0) = W(h = 1) = 0, which 
is impossible. It would be quite difficult to understand this result (the non-existence of a 
BPS bound and of solitons in the gauged model) without the additional insight provided by 
the superpotential equation. 

The case of the new baby Skyrme potential where the two boundary conditions W(h = 
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1.0 r 



0.8 - 




FIG. 3: Solution of the superpotential equation fll06j) for the new baby Skyrme potential 
V n = h(l — h), for the coupling constants A = 1, g = 1 and \i = 1 (W . . . dashed line). The 
solution develops a singularity at h ~ 0.9 and cannot be extended to the whole interval 

h G [0,1]. 



0) = and W(h = 1) = cannot be satisfied simultaneously corresponds to the generic case 
of two-vacua potentials. Still, there will exist some fine-tuned potentials for which the two 
boundary conditions can be satisfied (we shall see an explicit example in Section VI). These 
are precisely the fine-tuned cases where the extremum of V - which must exist because V has 
two vacua - coincides with the extremum of W - which must exist because W interpolates 
between W(h = 0) = and W(h = 1) = 0. It follows from the general results of subsection 
IV. C. 2 that also in these cases soliton solutions do not exist. We conclude that BPS solitons 
cannot exist for potentials with two vacua. 



6. The magnetic flux 

Here we want to demonstrate that the magnetic flux of a spherically symmetric soliton 
can be expressed in terms of the superpotential. Therefore, if we know the superpotential 
then we know the magnetic flux exactly. In the course of the derivation we shall also find the 
inequality a(y) > — 1 V y. The magnetic flux may be expressed in terms of the asymptotic 
value of a like 

$ = J rdrdipB = 27m J dya y = 27ma(yo) = 2ima 00 (H3) 

where yo is finite for compactons and infinite for non-compact solitons; further, n is the 
winding number of the spherically symmetric ansatz. Dividing the second BPS equation 
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(HOED by the first ffTUTD we find 



or 



where 



d y \n(l + a) = g 2 X 2 d y F (115) 

F ^w h - F <*>= 4 jf ,w K (116) 

which leads to the y integral 

\nC(l + a)=g 2 X 2 F(h(y)) (117) 

where C is an integration constant. Here we assumed that the potential is generic, i.e., that 
Wh = does not occur in the interval < h < 1, which implies that is finite in the 
same interval. At the vacuum h = 0, where Wh = 0, we assume that the potential behaves 
algebraically, i.e. V ~ h 2a for some a > 0, then Wh ~ h a , W ~ h a+1 near h = and 
is, in fact, zero at h = 0. As a consequence, F exists and is finite in the whole interval 
h G [0, 1]. For these generic potentials, it follows from the above result that 

a(y)>-l V y (118) 

and that the limit a — > — 1 may be reached only in the limit gX — > oo. The integration 
constant may be determined from the boundary conditions h(y = 0) = 1, a(y = 0) = 0, 

F(l) = -LlnC c = e 92x2pil) (119) 

9 X 

which, together with h(yo) = and F(h = 0) = leads to the asymptotic expression 

aoo = -i + e -9 2 ^ni) (120) 

which may be inserted into the expression for the magnetic flux. Specifically, in the limits 
of small and large electromagnetic coupling g we find for the magnetic flux 

g small : $ ~ -2nng 2 X 2 F(l) (121) 

g large : $ ~ -27m. (122) 

Both the small and large q behaviour coincide with the numerical findings for the full gauged 
baby Skyrme model in [25|, but here it is an exact result. We remark that this result is 
completely confirmed by our numerical calculations, as well. 
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V. NUMERICAL SOLUTIONS 

In this section we present the results of numerically solving the static field equations 
within the spherically symmetric ansatz ( 1361) . ( |371) . Concretely, we will consider the two 
cases of the old baby Skyrme potential V = 2h and of the potential V 2 = 2V S = Ah 2 . 

A. The old baby Skyrme potential 

The static field equations for the old baby Skyrme potential with the spherically sym- 
metric ansatz allow for a first integral which leads to the following equations (see Section 
2.4; y = r 2 /2) 

h y (l + a) 2 = ^(y-y ), (123) 

/ 2 \ 2 

(l + a)X y = \ 2 g 2 (-^ \ (y-y ) 2 , (124) 

where y$ is the integration constant of the first integral. Physically, yo is interpreted as the 
(squared) compacton radius. Further, the energy density is 

e = 2AV(1 + afh 2 y + 2fi 2 h + ^n 2 a 2 y , (125) 

and, taking into account the change of variable, dy = rdr, the static energy is (we choose 
the energy scale E — 1) 

E = J dye(y). (126) 

Obviously, the static field equations depend on the topological charge n and on the coupling 
constants A and /i only via the combination (fi/nX), therefore we may fix two of them and 
just vary the third one. Concretely we will choose n = 1 and A = 1, where the second 
choice just fixes our length scale. With this choice, // is a dimensionless parameter (coupling 
constant), and different values of \i correspond to different theories. 

Before performing the numerical calculations, it is useful to do a power series expansion 
about the compacton boundary. Inserting the expansion into the static field equations we 
find in leading order 

h = 4w2A2( f +6D)2 (y " vo) 2 + 0{g\y - y )% (127) 

a = b 0+ 12w4 / 2 2 ( f +6D)3 (y - Vo)' + 0{g\y - y Q f). (128) 

Here, bo is a free parameter (i.e., it remains undetermined by the equations). If we perform 
a shooting from the boundary, we have therefore the two free parameters yo and bo at our 
disposal which we may vary in order to satisfy the two remaining boundary conditions 
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Vo 


bo 


E 


Figure 


0.001 


6.325 


-1.1361 • 10" 6 


5.2996 


<m 


0.01 


6.324 


-0.134 • 10" 3 


5.2983 


0D 


0.1 


6.268 


-0.0135 


5.2794 


a 


1 


2.397 


-0.838 


3.6760 


® 


2 


0.810 


-0.9999964 


1.9711 


m 



TABLE I: Solutions of field equations for /i 2 = 0.1 and low g. 



h(y = 0) = 1 and a(y = 0) = 0. We have performed the shooting from the compacton 
boundary for the values fi 2 (/x 2 = 0.1, 1, 10, 50) with very similar results, so here we will only 
show the first one, /i 2 = 0.1. The main reason for this choice fi 2 = 0.1 is that this is the 
value chosen in reference 25j . We have also chosen different values of the coupling constant 
g. Thus, in Table [I] we show the parameter values and energies for the solutions found for 
the values (g = 0.001, 0.01, 0.1, 1, 2). Further, we show the graphs of the solutions in Figures 
H]- [6] for the values g = 0.1, 1,2. We do not display the figures for g = 0.001,0.01 because 
they look exactly like Fig. H] for g = 0.1, with the only difference that the graph for the 
gauge field a (and its derivative a') has to be multiplied by 10~ 2 for g = 0.01 and by 10~ 4 
for g = 0.001, because a is proportional to g 2 for small g to a high precision. 



On the other hand, for large values of g shooting from the boundary is problematic for the 
following reason. We cannot start the shooting exactly at the compacton boundary, because 
the fields take exactly their vacuum values at the boundary, and the numerical integration 
would only find the trivial vacuum solution and not the soliton. Instead, we have to start the 
integration slightly inside the compacton radius and use the above power series expansion 
for the determination of the "initial" value (i.e., boundary value). The problem is that this 
power series expansion in y — yo is, at the same time, a power series expansion in g. It is 
therefore reliable for small g but not for large g. The way out is to perform a shooting from 
the center for large g. A power series expansion in y about the center y = is again, at 
the same time, a power series expansion in g. The difference is that the fields do not take 
their vacuum values at the center, therefore we may start the shooting exactly at the center. 
Inserting the power series expansion at the center into the field equations, we get 

<y)-b iy + 9 ^y 2 + ... (130) 

so we have the free parameters b\ and yo to satisfy the two boundary conditions at the 
compacton boundary. We show the solutions for g = 5,10 in Table HT1 and in Figures [7] and 

m 



31 



9 


yo 


h 
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Figure 


5 


0.317 


-3.16261 


0.793768 


ra 


10 


0.158 


-6.33307 


0.397260 


m 



TABLE II: Solutions of field equations for /i 2 = 0.1 and high g. 



Now let us briefly comment on the behaviour of the soliton energies as a function of 
the gauge coupling constant g. We find numerically that the energies behave as 1/g for 
g > 1 whereas they are more or less constant when g < 1, for all values of fi 2 studied. 
The behaviour for small g exactly reproduces the analytical result found in Section IV. C, 
Remark 3. On the other hand, we were not able to find an analytic expression for the large 
g behaviour to compare with. The behaviour for different values of g for the constant value 
fi 2 = 0.1 is presented in Figure El In addition, studying the variation of the energy with \x 
for a fixed value of g, we find that it is proportional to /i (see Figure [TUl where g = 0.1). 

Another quantity of considerable physical interest is the magnetic flux 



rdrdtpB = 2nn / dya y = 27ma(y ) = 27m& 



;i3i) 



where in the case of the shooting from the boundary the magnetic flux may be expressed 
directly in terms of the free integration constant bo- We may infer from Table [I] that the 
magnetic flux grows like g 2 to a high precision for small g, whereas it approaches the constant, 
"quantized" value —2irn for large g. That is to say, the behaviour we find for the magnetic 
flux reproduces exactly the analytical results of Section VI. C. 6 and coincides with the one 
found in 25] for the full gauged baby Skyrme model. The behaviour of the magnetic field 
itself, on the other hand, is different from the results of 25j. We find that, for large g, 
a(y) changes almost exactly linearly from a(0) = to a(y ) ~ —1. The magnetic field 
B = na r /r = na y is, therefore, almost constant inside the compacton and rapidly decreases 
to zero near the compacton boundary y = yo (we remark that the (almost constant) value 
of a y for large g is too large to fit into the figures [7J El therefore we rescaled it by 1/10). 
Further, the compacton radius squared yo shrinks like g~ l for large g, so the compacton 
radius r = y/2y shrinks like g -1 / 2 , and the (constant) magnetic field inside the compacton 
grows like g. This should be contrasted with the findings of {25)], where they find a magnetic 
field which is almost completely concentrated in a tiny region about y = 0, i.e., r = for 
large g. 



B. Integrating the BPS equations 



In a next step, we want to compare the energies and solutions of the static field equations 
(I123p . (I124p with the energies and solutions of the BPS and superpotential equations, in 
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E B 


E 


0.001 


5.2984 


5.2996 


0.01 


5.2982 


5.2983 


0.1 


5.2783 


5.2794 


1 


3.6744 


3.6760 


2 


1.9705 


1.9711 


5 


0.793765 


0.793768 


10 


0.397259 


0.397260 



TABLE III: BPS bound and static soliton energies for /i 2 = 0.1. 

order to confirm that the solitons are, indeed, BPS solutions. In a first step, we compare the 
soliton energies with the BPS bound, where for the BPS bound we just have to determine 
W(h = 1), i.e., W(4>3 = —1) numerically for the old baby Skyrme potential, see eq. fl88|) . 
We have done this for different values of /z 2 with similar results, although, as before, here we 
only show the case /i 2 = 0.1. Then, in Table [TTT1 we show the values of the BPS energy, Eb, 
for each value of g, comparing them to the energies of our compacton solutions. In addition, 
in Figure [11] both energies are presented. 

The values of the BPS energies and the energies of our solutions agree with a precision 
of better than 5 • 10~ 4 in all cases, and all numerical energies are slightly above the BPS 
energies, as obviously must be true. We also remark that the values of Eb obtained for 
small g follow exactly the analytical expression. 

As a last step, we want to solve directly the two BPS equations together with the super- 
field equation, both as a consistency check of our numerical calculation and as a demonstra- 
tion of the BPS nature of the soliton solutions. We transform the superfield equation into 
a first order equation in the base space variable y by multiplying it by h 2 and by using the 
chain rule W y = Whh y and arrive at the system of three first order equations 

nh 2 y (l + a) + ^W y = 0, (132) 

na y + g 2 X 2 W = 0, (133) 

(the two BPS equations), and 

\w 2 + g 2 \ 2 h 2 W 2 - ^-h 2 y h = (134) 

(the superfield equation in base space), where y = r 2 /2. The boundary conditions we have 
to impose at the compacton boundary are 



W(y ) = , % ) = 



(135) 
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whereas 

W y {y ) = 0, h y (yo) = 0, a y (y ) = (136) 

are then consequences of the above equations. We want to solve these equations via shooting 
from the compacton boundary. Doing the power series expansion of the functions about the 
boundary y and imposing the corresponding boundary conditions, we find 

h = ^WTW iy - yo)2 + --- (137) 

■' yV •(.'/--.'/o!' • - (138) 



12n 4 A 2 (l + 6 ) 3 



We conclude that the expansions for h and a are exactly the same as above. This demon- 
strates that, in this case, we will indeed get exactly the same soliton solutions as above. 
The two free constants bo and yo are again used to implement the two remaining boundary 
conditions h(0) = 1 and a(0) = at the center. As a final check, we solve the system of 
equations ( 11321) . (I133p and (11341) numerically via shooting from the boundary. In Figure ( fl~2l) 
we present the solution for the superpotential for the case /i 2 = 0.1 with g — 0.1 (the graphs 
of h, a and the energy density are exactly like in Figure H]). The values of the constants we 
get for this solutions are basically the same as before. We find a similar situation for other 
values of fi 2 and (low) g. As explained before, we cannot use shooting from the boundary 
for high g because of the problems with the expansion at the boundary. 

C. The potential V = Ah? 

Here we want to briefly describe the numerical solution for the potential 

V = (1-n- 4>) 2 = Ah 2 . (140) 

The system of static field equations for the spherically symmetric ansatz is 

2u 2 

d y [h y (l + a) 2 ] - -£- 2 h = 0, (141) 
a yy =AX 2 g 2 (l + a)h 2 y (142) 

and the energy density reads 

e = 2AV(1 + a) 2 h 2 y + Afi 2 h 2 + ^n 2 a 2 y . (143) 

In this case the potential is quartic in small fluctuations about the vacuum, therefore we 
expect exponential-type solutions, so we have to shoot from the center. Performing an 
expansion about the center and imposing 

h(0) = 1, a(0) = 0, (144) 
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Figure 


0.1 


- 0.44498 


- 0.00892 


5.60583 


- 0.009940 


EH 


1 


-0.27393 


-0.707 


4.43440 


- 0.696595 


M 


2 


- 0.117166 


- 1.72637 


2.71178 


- 0.99763 


m 



TABLE IV: Solutions of field equations for ji 2 = 0.1 and g = 0.1, 1, 2. 

we get 

%) ~ 1 + f lV + - ^ y 2 , (145) 

a(y)^b iy + 2g 2 X 2 f 2 y 2 . (146) 

Now we have the two free parameters fi and &i to satisfy the two boundary conditions 
in the limit y — > oo. For the numerical calculation we choose n = 1, A = 1, /i 2 = 0.1 
and g = 0.1, 1,2. Then we find soliton solutions for the values shown in table HVl On the 
other hand, solving the superpotential equation (which has a global solution in this case) 
numerically, we find the BPS bound Eb = 2ir\W(h = 1)| which agrees with the soliton 
energies E within the shown precision in all three cases, therefore we do not show them 
separately. Apparently, the numerical convergence is even better for this potential. The 
result of the numerical integration of the soliton solutions is shown in Figures [13], [HJ [151 
For g = 0.1 and g = 1 the exponential approach to the vacuum is clearly visible. For g = 2, 
on the other hand, the approach to the vacuum seems to be more like a compacton, and one 
wonders whether the exponential approach continues to hold for large g. We shall find in the 
next Section that the approach is, in fact, exponential, and the compacton-like appearance 
is due to a very fast exponential decay for large g, essentially like exp(— e g2x2 y). 

Finally, we briefly comment on the magnetic flux and magnetic field for this case. Like 
in the case of the old baby Skyrme potential (Section V.A) and in Ref. {25}, the magnetic 
flux $ = 27ma OQ grows like g 2 in absolute value for small g and approaches the constant 
value —2im in the limit of large g, in accordance with the exact results of Section IV. C. 6. 
Further, the magnetic field itself in the limit of large g behaves like the one in Section V.A, 
i.e., it is almost constant in the core of the soliton and drops to zero quickly in a thin shell 
of rapid exponential decay, see Fig. [T5j 



VI. EXACT SOLUTIONS 

Here we want to show that for some fine-tuned choices of potentials there exist exact 
solutions both for the superpotential equation and for the corresponding BPS equations. It 
is in fact easy to find exact solutions to the superpotential equation. We just start with 
a given superpotential and calculate the potential from the superpotential equation (1106!) . 
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The disadvantage of this method is that it somewhat obscures the physical meaning of the 
coupling constants, because the potential now depends on them. For some fixed, given 
values of these constants, on the other hand, this is a perfectly valid method to produce 
exact solutions. Let us choose as a first example 

W = -U 6 (147) 
X z 

where b > 1, then the potential is 

nV^^h^^ + fXVy (148) 

Besides being interesting in its own respect, this exact solution may help us to better un- 
derstand the numerical solutions for potentials of the type V ~ h 2a studied in the previous 
Section, for the following reason. For small g, the first, g independent part of the potential 
(I148P dominates and, therefore, the solution W of the superpotential equation for the poten- 
tial V ~ h 2a may be approximated by the exact solution ( I147p . If we are interested only in 
the asymptotic behaviour, then the restriction to small g can be lifted, because sufficiently 
close to the vacuum (i.e., for sufficiently small h), the first term (proportional to b 2 ) in the 
potential (I148p dominates even for large g. We shall see below that for the potential (I148P 
not only the superpotential is known, but also the spherically symmetric BPS equations 
can be solved exactly. We can, therefore, use the asymptotics of these exact solutions to 
determine exactly the asymptotics of solitons for the potentials V ~ h 2a , and especially for 
the potential of Section V.C. 

Before doing so, let us briefly mention that by this method we may also find global 
solutions for the superpotential equation for two- vacua potentials. Indeed, choosing, e.g., 
the superpotential 

W = -U 2 (l - h) 2 (149) 
X 2 

we find the potential 

V = ^ 2 h 2 (l - h) 2 [(1 - 2h) 2 + g 2 X 2 h 2 (l - h) 2 ] (150) 

where, as explained in Section IV. C. 2, the extrema of V and W coincide. It follows from the 
results of Section IV. C. 4 that, still, BPS soliton solutions do not exist for this potential. Be- 
sides, the existence of a global solution W for this two-vacua potential V is the consequence 
of a fine-tuning of the parameters of V. 

A. Exact BPS solutions 

We shall find several examples where for an exact expression for the superpotential the 
spherically symmetric BPS equations 

2nhy{l + a) = ~W h (151) 
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na y = -g 2 X 2 W (152) 

have exact solutions, too. 

1. A first example 

As a first example, we choose the superpotential W = -kh 2 , with BPS equations 



1 



2nh y (l + a) = -— h (153) 
A 

na y = -g 2 h 2 . (154) 



Introducing p = h 2 we get 

p_ 

A 2 

na y = —g 2 p. (156) 



np y {\ + a) = -— (155) 



Inserting the second formula into the first one gives 



and from (11561) 



Then, 
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Pw = —^ a vv ( 158 ) 



no ro = -^d y ln(l + a) =>- -nd yy (l + a) = --r<9 a ln(l + a) (159) 
A A 



which leads to the first order equation 



d y (l + a) = ^ln[C(l + a)]. (160) 



It can be further integrated 

1 f d(C(l + a) 1 



CJ ln[C(l + a)] n\ 2 
where B and C are integration constants. Thus finally 



(y - B) (161) 



Lu[C{l + a{y))] = ~{y-B) (162) 

where Li is the logarithmic integral function. The general solution for the profile function h 
may be derived from (11561) and (1 1 6 [) 

h 2 = ^- 2 \n[C(l + a)}. (163) 
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It remains to determine the integration constants from the boundary conditions for a soliton 
solution. From h(0) — 1, a(0) = 0, and Eq. (I163p we find 

_ In (7 

whereas from a(0) = and Eq. (11621) we get 



C = e 9 ^ (164) 



Li [C(l + a(y))} = Li [1] = -— -(y - B). (167) 



^-UC=^- => B = nAV^Li (e 9 ^ 2 ). (165) 
C n\ 2 

In order to find whether the fields take their vacuum values at a finite or infinite yo, we first 
insert the boundary condition h(y ) = into Eq. (1 163[) to conclude 

C7(l + a(y )) = l (166) 
and then use this result in Eq. ( 1 162ft to find 

C 
^X 2 

It is one of the properties of the logarithmic integral that Li [1] = — oo, therefore we conclude 
that yo = oo. We may also determine the asymptotic value of a(y), 

a(y = oo) = = -1 + CT 1 = -1 + e~ g2x2 (168) 

which implies doc ~ — g 2 X 2 + 0(g 4 X A ) for small g and ~ —1 for large g, exactly repro- 
ducing the numerical findings of Section V.C (we remind the reader that the asymptotic 
behaviours of the exact solution of this section and of the numerical one of Section V.C are 
the same). 

To summarize, the soliton solution reads 

<y) = - 1 + ^ u ~ 1 ( uc -§?y) ( 169 > 



(170) 



with C given in (j!64p . The approach to the vacuum may be determined from the asymptotic 
behaviour of the logarithmic integral near one, 

Li(x) fa In |1 — x\, x — > 1. 

Hence, at y — > oo we have that (7(1 + a(y)) — > 1 and from (j!62p 

L\n\l-C(l + a(y))\ = -^, y -> oo (171) 

Thus, 

o(y) ^ -1 + - + - e ~^ = aoo + e~ 9 exp -f , (172) 



C C 00 r I nA 2 
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(i73) 

where the last result comes from the second BPS equation na y = —g 2 h 2 . It follows that these 
solitons (as well as the ones of Section V.C) are exponentially localized, as was announced 
already in Section V.C. 



2. More examples of BPS solutions 

The BPS equations can, in fact, be solved for all potentials f)148p . i.e., for all superpo- 
tentials ()147p . The BPS equations are 

2nh y (l + a) = T ^h b - 1 (174) 

na y = T g 2 h b (175) 

where, depending on the particular case, one has to take the plus or minus sign. Repeating 
the previous computations we get 

h 2 = b ~HC(l + a)] (176) 

^i(j^y\ y _ B) (m) 



C J (ln[C(l + a)]) b / 2 ^ n \2g 2 X 
The boundary conditions at y = give 

\nC = ^X 2 g 2 . (178) 

The boundary condition at y = yo gives again the asymptotic value of the magnetic function 
a(yo) 

C(l + a(y )) = l a(yo) = e-^ 2 -l (179) 

In order to calculate yo and B one has to compute the integral, which strongly depends on 
the parameter b. For example for b = 4 we get (we take the plus sign) 

The l.h.s. function Li(x) — is a function which starts at zero value at the origin and then 
grows to infinity at x — 1. Hence, yo = oo i.e., it is a non-compacton configuration. Finally, 
as a(0) = we get 

C 1 / 2 x 1 



^C- — = --[ — ) B (181) 



lnC n \gX 

or 
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VII. CONCLUSIONS 

In this article we investigated in depth the gauged version of the BPS baby Skyrme model, 
i.e., the BPS baby Skyrme-Maxwell system. We found that, like the ungauged model, the 
gauged model still has infinitely many symmetries, both Noether and non-Noether ones. 
Further, it continues to have a BPS bound and to support BPS soliton solutions saturating 
this bound. The BPS bound is, however, quite different from other, known BPS bounds. 
Known BPS bounds either bound the energy (or euclidean action) in terms of the topological 
degree only (e.g., instantons, 0(3) nonlinear sigma model) or involve both the topological 
charge density and the potential of the theory (e.g., scalar field domain walls, or the ungauged 
BPS baby Skyrme model). Here, on the other hand, we find a BPS bound in terms of the 
topological charge density and an auxiliary function which is related to the potential via 
a first order differential equation. We called the auxiliary function " superpotential" and 
its defining differential equation "superpotential equation" due to their similarity with the 
corresponding equation in supergravity, as explained Section IV.A. A further consequence 
of this new type of BPS bound is that it is a rather nontrivial problem whether for a given 
potential BPS soliton solutions exist, at all. One necessary condition is the existence of the 
superpotential on the whole target space (i.e., a global solution of the superpotential equation 
( I106jl ). but we found examples where BPS soliton solutions do not exist despite a globally 
existing superpotential, so this condition is, in general, not sufficient. Specifically, we found 
that for potentials with more than one vacuum BPS solitons never exist. The ungauged 
model with multi-vacua potentials, on the other hand, supports both a BPS bound and 
BPS solitons, so we find the interesting result that for systems where the matter sector 
does not have a unique vacuum, the inclusion of the electromagnetic interaction destabilizes 
the "particles" (topological solitons) of the theory. For models with monotonously growing 
potentials, instead, we found both a nontrivial BPS bound and BPS soliton solutions in all 
cases we considered, motivating the conjecture that this will always be true. 

We studied rotationally symmetric solitons numerically for some concrete potentials and 
found in all cases that the soliton energies, indeed, saturate the corresponding BPS bound. 
We also found that the soliton energies decrease with increasing electromagnetic coupling 
(e.g., like g -1 for the old baby Skyrme potential, see Section V.A). If the superpotential is 
known exactly, we found that, at least in some cases, the whole system of BPS equations can 
then be solved exactly, providing us with analytic, explicit expressions for the solitons, see 
Section VI. A further quantity which can be calculated exactly once the superpotential is 
known is the magnetic flux, see Section IV. C. 6. Finally, we want to point out that the BPS 
bounds we found for the gauged BPS baby Skyrme models for different potentials provide, 
at the same time, BPS bounds for the corresponding full gauged baby Skyrme models which 
will, in many occasions, provide tighter bounds that the bounds already known (see our 
discussion in Section II. B). 

The model has shown a surprisingly rich mathematical structure and there remain many 
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interesting open problems which deserve further investigation. One first, obvious issue is to 
better understand under which conditions BPS soliton solutions will exist. That is to say, 
to find rigorous mathematical answers to the following two questions, 

i) for which potentials V(h) the superpotential equation fll06[) has global solutions in the 
whole interval h 6 [0, 1] (i.e., on the whole target space), and 

it) which additional conditions the potentials have to obey such that BPS soliton solutions 
exist. 

Another question of interest is whether the BPS property of the model may be related 
to some further structure as, e.g., supersymmetry, as is frequently the case. In this respect, 
already the ungauged model involves a surprise. Indeed, while it was found recently that 
the full baby Skyrme model allows for a supersymmetric extension 53j, this supersymmetric 
extension is not possible for the restricted (i.e., BPS) baby Skyrme model, which indicates 
that the relation between BPS equations and supersymmetry is more involved for field 
theories with non-standard kinetic terms. 

From a more physical point of view, the most important question is, of course, if and 
to which degree the structure found in this model can be generalized to the Skyrme model 
in 3+1 dimensions. A first important observation is that, as was already discussed in the 
introduction, the ungauged Skyrme model, too, has a submodel (the BPS Skyrme model) 
sharing all the nontrivial features with the BPS baby Skyrme model (integrability, BPS 
bound and exact BPS solutions). In this sense, the hope to be able to generalize some of 
the results of the present paper to 3+1 dimensions is well-founded. The gauged version of 
the BPS Skyrme model will, nevertheless, most likely present additional difficulties, where 
the most obvious one is related to the fact that the magnetic field is a pseudoscalar in 
2+1 dimensions, whereas it is a pseudovector in 3+1 dimensions. A further question of 
interest will be whether the gauged model in 3+1 dimensions maintains all the symmetries 
of the ungauged model, as is the case in 2+1 dimensions. In any case, this problem shall be 
investigated in future publications, where we hope that the detailed analytical and numerical 
results developed in the present article should be instrumental in the understanding and 
investigation of the 3+1 dimensional system. 

Finally, another question of considerable interest is, in our opinion, whether BPS bounds 
of the type found in the present article can be applied in a more general context, i.e., outside 
the field of Skyrme type models. 

Note added: While finishing this paper we became aware of the preprint (54)], where 
the BPS bound we derived in Section IV.A was found using completely different methods, 
extending previous results of the so-called Bogomolny decomposition for the ungauged model 



55|. 
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FIG. 4: Solutions for g = 0.1. 
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FIG. 5: Solutions for g — 1. 




FIG. 6: Solutions for g = 2. 
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(a) Function h (magnified near the compacton boundary) and 
its derivative (dashed line). 
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(b) Function a and its derivative (rescaled by 10 1 ; dashed line). 
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FIG. 7: Solutions for g = 5. 
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(a) Function h (magnified near the compacton boundary) and 
its derivative (dashed line). 




(b) Function a and its derivative (rescaled by 10 1 ; dashed line). 
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FIG. 8: Solutions for g = 10. 
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FIG. 9: Static energy as a function of g for fi 2 = 0.1. The continuous line shows the 
behaviour 1/g for g > 1. We have found a similar behaviour for the other values of /x 2 . 
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FIG. 10: Static energy as a function of \x for g = 0.1. The continuous line shows the 
behaviour proportional to fi. We have found similar plots for the other values of g. 
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FIG. 12: Solutions for the superpotential of the BPS system for /x 2 = 0.1 and g = 0.1 

(W . . .dashed line). 
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FIG. 13: Solutions for fi 2 = 0.1 and g = 0.1. 
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FIG. 14: Solutions for jj 2 = 0.1 and g = 1. 
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